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Using the perturbation method based on a variational phonon basis obtained by the modified 
Lang-Firsov (MLF) transformation, the two-site single polaron Holstein model is studied in pres- 
ence of a difference in bare site energies (e<j=e2-ei). The polaronic ground-state wave function is 
calculated up to the fifth order of perturbation. The effect of ed (acts as a site-energy disorder) 
on the polaron crossover, polaronic kinetic energy, oscillator wavefuncion and polaron localization 
are studied. Considering a double-exchange Holstein model with finite e^, role of disorder on the 
properties of the double-exchange system is also discussed. 
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I. Introduction 

Study of narrow band electronic systems with strong 
electron-phonon (e-ph) interaction has long been an ac- 
tive research area in condensed matter physics. The field 
has drawn renewed interest following evidence of pola- 
ronic charge carriers in underdoped high-T c cuprates , 
manganites 2], and organic superconductors. The one- 
dimensional polaron problem is also relevant in semicon- 
ductor physics, quantum dots and linear conjugated 
organic polymer conductors p|. The simplest model for 
studying polarons is the Holstein model where an 
electron in a narrow band interacts locally with optical 
phonons. For large e-ph coupling the polaron is a small 
polaron with high effective mass, while for small coupling 
it becomes a large polaron having a much lower effec- 
tive mass for a finite adiabatic parameter. The crossover 
from a large to a small polaron and the corresponding 
change in the polaronic properties in the ground state 
have been studied for the Holstein model by different 
groups Q) IE 13 using various methods enlightening 
our understanding in this field. However studies on the 
nature and properties of polarons in presence of disorder 
are few and need much more attention. The imperfec- 
tions or disorder may play an important role in complex 
materials (high-T c oxides, manganites etc.) where signa- 
tures for polaronic carriers are found. Recently the small 
polaron concept has been used to explain the charge mo- 
tion in DNA where the electronic band is very narrow and 
the presence of different kind of molecular units induces 
large disorder potential 10]. 

In absence of any disorder, translational symmetry en- 
sures that the polaronic ground state is delocalized how- 
ever large is the e-ph coupling strength provided other 
parameters (electronic hopping, phonon frequency) are 
finite. The large to small polaron crossover is a contin- 
uous one which is consistent with the ground state 
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properties being analytic functions of e-ph coupling 
Localization requires a breakdown of the translational in- 
variancc which may be achieved through randomness of 
the site potential or hopping. The effect of site diago- 
nal disorder on polaronic properties has been addressed 
by some authors. Shinozuka and Toyozawa |12| studied 
disorder induced self-trapping in a tight binding model 
in which the local site energies are randomly distributed 
between two values and found that the exciton-lattice 
interaction acts with the disorder to produce severe lo- 
calization associated with a self-trapped exciton. In his 
study the lattice vibration was treated as classical os- 
cillators. Bronold et. al |13| studied similar model but 
with an infinite coordination number within the dynam- 
ical coherent potential approximation. However the lim- 
itation of the coherent potential approximation is that it 
cannot fully distinguish between localized and itinerant 
states. Bronold and Fehske 0] improved the method to 
overcome the above shortcoming. They followed statis- 
tical dynamical mean field theory to predict localization 
of small polarons by extremely small disorder. However 
a proper study of the effect of disorder on the polaron 
crossover is not made in any of the above investigations. 
In the present work we consider a two-site cluster with 
different site energies. This is the minimal system to 
study the competition of the inter-site electronic hop- 
ping with the localization induced by the combined ef- 
fect of the e-ph coupling and the site energy disorder. 
Difference in site energies would remove the two-fold de- 
generacy of the system in the absence of hopping and 
would tend to localize the electron in the lower potential 
site. For convenience we will refer to the difference in 
site energies as disorder strength because it partly mim- 
ics the role of disorder in larger systems. Another aspect 
of choosing a two-site Holstein model is that almost ex- 
act results may be obtained for such a system with the 
perturbation method [l5| using a modified Lang-Firsov 
(MLF) basis [llj. For the Holstein model the interaction 
is very short-ranged and the essential physics relating po- 
laronic behavior for a larger system is similar to that in a 
two-site system. In studying a Hubbard-Holstein model 
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similar conclusion has been reached in Ref. 17] . The 
relevance of studying a two-site system in the context of 
Holstein and Holstein-double exchange models has been 
discussed in details in Ref. [Isj . 

In Sec. II we discuss the formalism and perturbation 
calculations. In Sec. Ill we present the results obtained 
by MLF method and discuss the role of the disorder 
strength on the polaron crossover and the kinetic energy 
of the system. The localization of the polaron and ground 
state polaronic wavefunction will also be discussed. Ex- 
tension of this model to the double-exchange model is 
included in Sec. IV. Finally, a summary of the results 
are presented in Sec. V. 



II. Formalism 

The two-site single-polaron Hamiltonian is 

H = emi + e 2 n2 - ^2t(c\ a C2a + c' 2a c 1(T ) 

a 

+ gw^m^bi + b^ + u^blbi (1) 

i,a i 

where i —I or 2, denotes the site. t\ and e 2 are the 
bare site-energies at site 1 and 2, respectively. Cj ff (c\ a ) 
is the annihilation (creation) operator for the electron 
with spin a at site i and rii a (=<J irT Ci a ) is the correspond- 
ing number operator, g denotes the on-site e-ph coupling 
strength, t is the usual hopping integral, bi and b\ are 
the annihilation and creation operators, respectively, for 
the phonons corresponding to interatomic vibrations at 
site i and to is the phonon frequency. This Hamiltonian 
has spin degeneracy for the one electron case so the spin 
index is redundant. 

Introducing new phonon operators a = (b\ + fe 2 )/-\/2 
and d = (b\ — 6 2 )/-\/2, the Hamiltonian is separated into 
two parts (H = Hd + H a ) : 

Hd = emi + e 2 ri2 - t(c\c2 + c|ci) 

-I- uig+{ni - n 2 ){d + S) +wcftd (2) 
and H a = ua'a — ujn 2 g+ (3) 

where g+ = g/V%, a = a + ng + . The a-oscillator couples 
only with the total number of electrons n(— n\ + 712), 
which is a constant of motion. Consequently H a de- 
scribes just a shifted oscillator, while Hd represents an 
effective e-ph system where phonons couple with the elec- 
tronic degrees of freedom. For a perturbation method it 
is desirable to use a basis where the major part of the 
Hamiltonian becomes diagonal. When the hopping (t) 
is appreciable compared to the polaron energy (g 2 uj), a 
retardation between the electron and associated lattice 
distortion sets in. The retardation is important for low 
and intermediate values of g; it produces a spread in the 
polaron size, the resultant polaron is a large one. In Ref. 



[IH we have shown that the MLF perturbation method 
works much better than the Lang-Firsov (LF) method 
for a large region of parameter space where the retar- 
dation is important. In the strong-coupling limit the 
MLF method reduces to the LF method and it works 
well there. We use the MLF transformation where the 
lattice deformations produced by the electron are treated 
as variational parameters 0, . For the present sys- 
tem, Hd = e R Hde~ R where R = X(ni — n 2 )(rf t — d) and 
A is a variational parameter describing the displacement 
of the d oscillator. 

The transformed Hamiltonian is then obtained as 

Hd = Lod)d + (ei — e p )m + (e 2 — e p )n 2 

- t [c\c 2 exp(2A(dt - d)) 

+ c 2 ci exp (-2A(d f -d))] 

+ Lu(g + - X)(m - n 2 ){d + d ] ) (4) 

where e p — w(2g + — A)A. 

For the single polaron problem we choose the basis set, 

\+,N) = {a lC \+a2cl)\0) e \N) ph 

\-,N) = (a 2 c t 1 - ai c 2 )|0) e |7V) ph (5) 

where |+) and |— ) denote the electronic states and \N) 
denotes the iVth excited oscillator state in the MLF 
phonon basis. The normalization condition requires 
af + a 2 , = 1 • The unperturbed part of the Hamiltonian 
is chosen as 

Ho = (jd'd + (ei - e p )ni + (e 2 - e p )n 2 

- t e (c[C2 + cjci) (6) 

where t e = t exp(— 2 A 2 ). The remaining part Hi = (Hd — 
H ) is treated as a perturbation. The states |±,iV) are 
the eigenstates of the unperturbed Hamiltonian (Ho) for 

r = 01/02 = [fd + V td 2 + 4£ e 2 ]/2i e where e<i = e 2 - €\. 
The corresponding eigen energies are given by 

E% = Nu, + _ £p T 1 ^/ (ei _ C2 )2 +4te 2 (7) 

The state |+, 0) has the lowest unperturbed energy, 

4° } = + T " e P ye d 2 +At 2 . 

The general off-diagonal matrix elements between the 
states \±,N) and \±,M) are calculated for (N — M) > 
as follow : 
for even (N - M), 

(JV.il/Txli, M> = T-i e7T ^T (8) 
(1 + H) 

{N,±\H X \^M) = -te^^- (9) 
(1 + r A ) 

for odd (N-M), 

(jV,±|£ri|±,A0 = tVNu>( 9+ r l Sn,m+i (10) 
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{N,±\Hi.\T,M) = ±t e + VNuj(g+-X) 



2r 



(1 + r 2 ) 



(11) 



It may be noted that for the ordered case (ed—0, hence 
r=l) the off-diagonal matrix elements (N, ±|i/i|±, M) 
are nonzero only for even (N — M) while (N,±\Hi\^p, M) 
are nonzero only for odd (N — M) . 

Now one has to make a choice of A so that the pertur- 
bation corrections are small and perturbative expansion 
becomes convergent. Here we will follow the procedure of 
our previous work |l5j| to find out the variational phonon 
basis as a function of e-ph coupling. Minimizing the un- 
perturbed ground state energy E^ with respect to A we 
obtain 



A = 



(12) 



The perturbation corrections of different orders to the 
ground-state energy and the wavefunction may be calcu- 
lated using the general off-diagonal matrix elements in 
Eqs.(S-ll). The ground-state wave function in the MLF 
basis may be written as, 



\G) = 



1 



,0) 



N=l,2, 



-,N) 



E c n\ 



-,N) 



(13) 



The coefficients are obtained in terms of the off- 
diagonal matrix elements and unperturbed energies and 
Nq is the normalization factor. 

The correlation functions involving the charge and the 
lattice deformations (n 1 u 1 ) and (n. 1 u 2 )o> where U\ and 
U2 are the lattice deformations at sites 1 and 2 respec- 
tively, produced by an electron at site 1, are the standard 
measure of polaronic character and indicate the strength 
of polaron induced lattice deformations and their spread. 
Following Ref. the correlation functions may be writ- 
ten as 



(nxu 2 ) 
where Aq 



<9+ + X) Nc- + Nc- 



<9+ -A) 



Bo_ _Ao_ 

N G N G 



(G\m(d + (t)\Cf) = 



(14) 
(15) 

-(rcj~ + r 2 c+) 



(1 + r 2 ) L 



.2„+, 



+ r(c N c 



+ 1 "T" ^N^N+l) + c N c N+l\ 



and B = (G\ ni \G) 



? 7TT>, 



(1 +r 2 ) ^ (1 + r 2 ) 

v 1 N=l v ; 



\2rc N c+ +r 2 \c+\' 2 + \c, |2 



The static correlation functions in Eqs. (14) and (15) 
are calculated for different to examine the role of the 
disorder on the polaron crossover. 

The kinetic energy of the system in the ground state 
is obtained as 

Ek.e = <G\H t \G> 



1 

Ne~ 



2t e axa 2 + 2 E E c ^ e ' N \ H t\°> +> 

e=± N^Q 



E E c%c e M {e,N\H t \M,e') 



(16) 



where H t = —t[c\c2 exp(2A(c^ — d)) + c3,ci exp(— 2A 
(d^ — d))] is the kinetic energy operator in the MLF basis. 
The occupation number n(k) of the charge carrier for the 
ground state within MLF method is also calculated as 



i(4±4)( Cl ± C2 ) 



(G\noM = 



1± 



E 



K.E 



t 



(17) 

(18) 



where 0, it denotes the values of the wave vector k. There- 
fore, it basically reflects the nature of the kinetic energy. 

The ground-state wave function for the d-oscillator is 
obtained from Eq. (13) by using the wavefunction for 
the N-th excited (MLF-displaced) harmonic oscillator for 
\±,N). 



ip n (x) = (x\n) 



1 



7r 1 /4^2"n! 



where H n (x) is the Hcrmite polynomial of degree n and 
xq is the displacement of the oscillator due to the MLF 
transformation. 

To study the localization effect due to disorder in site 
energy we calculate | e ( 1 1 G) | 2 and | e (2 1 G) | 2 . These are the 
probabilities that the polaron in the ground state (|G)) 
lies at site 1 and site 2, respectively. 



^ - Wo 



K| 2 + E \ aiC N + a 2 c Af| 2 



|fl2| 2 + E \ a2 °N ~ a l C Jvl 5 
JV#0 



(20) 
(21) 



III. Results and discussions 



In this paper all the results are derived by calculating 
the ground-state wavefunction up to the fifth order of 
perturbation. For the ordered case a comparison of our 
MLF-perturbation results 0] with the exact results by 
Rongsheng et. al [2(j shows that the MLF method up to 
the fifth order gives exact results for t — 0.5 (in a scale 
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FIG. 1: Plot of the correlation function (n\U2) versus g+ for 
t/u = 0.5 for different values of e^. Inset: on-site correlation 
function (niui). 



where ui=l) whereas for higher values of t (=1-1 and 2.1) 
very accurate results are produced by the MLF perturba- 
tion method for both strong and weak coupling regions. 
However in a narrow region of intermediate coupling the 
perturbation results for high values of t deviate from the 
exact results. In this paper we present the results mainly 
for the nonadiabatic regime (t/u < 1.0) for which the 
convergence of the MLF-perturbation series for both or- 
dered and disordered cases is found to be very good in 
the entire region of the e-ph coupling strength. For the 
ordered case, the convergence in energy and correlation 
functions has already been reported in Ref. and that 
for the wave function in Ref. [2l| . 

In Fig. 1 we plot the variation of the correlation func- 
tions (riiUi) and {niUz) with g + for t/cu = 0.5 and for dif- 
ferent values of the disorder strength e^. For intermediate 
coupling {riiii?) has appreciable values and (niui) / (ni) 
(shown in inset) deviates downwards from its small po- 
laronic LF value of 2g + . These are the signatures of 
retardation and directly show the spread up of the po- 
laron. With increasing the value of (111U2) decreases 
and {n\U\) shows less deviation from its LF value. This 
points to a reduced retardation effect with increasing dis- 
order strength. 

The polaron crossover from a large to a small po- 
laron may be studied by the correlation function \ = 
(ni(iti — U2)) /2g + (rii) . In the small polaron limit the 
retardation effect is negligible and x S e ts its standard 
LF value (=1). For a larger size polaron the value of \ 
is lower. In Fig. 2(a) we plot x as a function of e-ph 
coupling strength. It is seen that the size of the polaron 
becomes more localized with increasing disorder strength. 
The polaron crossover is continuous, but a change in the 
curvature of x vs - 9+ plot is observed at a point in the 
crossover region. This point (in g + space) where the cur- 
vature in x vs - 5+ plot changes may be taken as the 
crossover point. The crossover point shifts to lower value 
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FIG. 2: (a) The variations of x — (ju(ui — 112)) /2g+(m) with 
g+ for different values of td for t/u) = 0.5. (b) The variation 
of the kinetic energy Ek.e with g+ for t/uj—0.5 and ed~0, 
0.5 and 1.0. (c) Different parts of the kinetic energy (see the 
three terms in Eq (16) referred as K.E(l), K.E(2) and K.E(3), 
respectively) vs. g+ for different e^. 



of g + as €d increases. Thus the disorder favors formation 
of small polarons. In Fig. 2(b) we plot the kinetic energy 
of the polaron as a function of g+ . The kinetic energy is 
suppressed in presence of disorder in the range from low 
to intermediate values of g+. But disorder has no signifi- 
cant effect on the kinetic energy for strong coupling where 
the polarons are small. The kinetic energy in Eq. (16) 
has three terms. The first term is the contribution to the 
kinetic energy from the unperturbed ground state which 
is significant for low to intermediate couplings. The sec- 
ond term is the contribution from the matrix elements 
of Ht connecting the unperturbed ground state to higher 
phonon states. The third term is due to the hybridiza- 
tion of higher phonon states through H t . We plot these 
contributions separately as a function of g + in Fig. 2(c) 
for £^=0 and 0.5. It is found that the first contribution, 
which is dominant in the range from low to intermediate 
couplings and proportional to the product of the coeffi- 
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MLF (E d =0) 
exact (£d=0) 
MLF ( ia=0.5) 



t = 0.5, 00=1.0 




FIG. 3: The wave vector dependence of the occupation num- 
ber n(k = 0, 7r) of the charge carrier for the ground state with 
e-ph coupling within MLF method for 6^=0 and 0.5. Exact 
results (for e<j=0) from Ref. |22T| is also plotted for compari- 
son. 



cients a± and ct2 (see Eq.(16)), is suppressed substantially 
with increasing disorder strength. For strong-coupling 
regime the contribution to the kinetic energy comes only 
from the second term of Eq.(16) and that one is not af- 
fected by disorder in this region. It may be mentioned 
that contribution from the third term is too small to show 
up in the figure for the range of our study. 

In Fig. 3, we have plotted the occupation number n(k) 
of the polaron for the ground state with e-ph coupling. 
The difference between n(fc=0) and n(k=ir) depends di- 
rectly on the delocalization of the polaron (see Eq (18)). 
In absence of hopping this difference vanishes. With in- 
creasing e-ph coupling the difference between n{k) re- 
duces as the kinetic energy is more and more suppressed. 
For the disordered case and in the range from weak to 
intermediate coupling, the difference in occupations re- 
duces further owing to disorder-induced suppression of 
the kinetic energy. In the strong-coupling regime, n(k) 
does not show any change with ea because the kinetic 
energy is almost independent of in this region. The 
exact results for the occupation number (2^ for the or- 
dered case (ed=0) are also shown in Fig. 3 for comparison 
with our MLF results. 

In Figs. 4(a-d) we have shown the ground-state wave 
functions of the d oscillator for different values of e-ph 
coupling for both the ordered and disordered cases con- 
sidering that the electron is located at site 1. The re- 
sults for the ordered case have been discussed in previ- 
ous works [ijj, 1 2(1 l22j . The main changes, which are 
observed for the disordered case (ed=1.0) in the antiadi- 
abatic regime, are as follow: 

(i) for weak coupling (g + = 0.4) the wave function, which 
shows displaced Gaussian-like single peak, is slightly 
more shifted; (ii) for intermediate coupling (g + = 1.3), 
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FIG. 4: Oscillator wave function (not normalized) vs. x with 
tjw = 0.5 for e d = and 1.0 for (a) g+ = 0.4, (b) g+ = 1.0, (c) 
g+ — 1.3 and (d) g+ = 2.0. (e) The coefficients (see Eq. 
(13)) of the ground state wave function with phonon number 
N for t/u>= 0.5 and g+=2 for ed=0 and 1.0. Solid curve gives 
the value of for odd N and for even N for the ordered 
case. For disordered case the thick dashed curve shows 
and the thin dashed curve shows c~^. (f) Magnitude of the 
coefficients (|cjv|) are shown for the same case as (e). 



the additional shoulder, which appears (at the right side 
of the main peak) for the ordered case, is absent in pres- 
ence of disorder; (iii) for strong coupling (g + = 2.0) 
though the wavefunction around the main peak is iden- 
tical for both cases, the additional small broad peak, ob- 
served for the ordered case, is not found for the disor- 
dered case. The first feature indicates that the on-site 
polaronic deformation is larger for the disordered case; 
while the second feature may be due to the reduced re- 
tardation effect. The third feature, i.e., disappearance 
of the additional broad peak for the disordered case, has 
a contradiction with the behavior of the kinetic energy. 
Both this broad additional peak and the kinetic energy in 
the strong-coupling regime are determined by the coeffi- 
cients c N in Eq.(13). A logical question then arises why 
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FIG. 5: Plots of | e <l|G)| 2 and |e(2|G}| 2 versus g+ for t/w = 
0.5. | e (l|G)| 2 and | e {2|G)| 2 are the probabilities of site 1 and 
site 2 being occupied by the electron, respectively, (a) The 
solid line is for the ordered case (e<j=0) where both occupancy 
are same, thick dashed line for ed=0.1 and thin dashed line for 
£d=0.5. (b) The same curves for £^=0.001. Polaron crossover 
(x) is also shown, (c) Plots of x an d |e(l|G)| 2 for different 
values of t. x is denoted by different symbols and | e (l|G)j 2 
by broken line. 



the broad peak in the oscillator wave function is modified 
by disorder while the kinetic energy in that region is un- 
affected. In Fig. 4(e) we plot versus N for g + = 2.0 
for both ordered and disordered cases. For the ordered 



case c N is nonzero for odd N, while cjy is nonzero for 



■N 



even N. For the disordered case both the coefficients 
and c^r are nonzero for any N. For strong coupling, the 
reduced polaronic hopping makes the value of r (see the 
expression of r after Eq.(6)) very large in case of disorder. 
As a consequence we find that takes negligible value 
except for N = 1, while cjj is appreciable for both odd 
and even AT. changes sign alternately with N . How- 
ever, the magnitude of the coefficients (|c^|) for any N 
do not change much with disorder strength in the strong- 
coupling regime (shown in Fig. 4(f)). In this region we 



observe that the kinetic energy essentially depends on the 
magnitude of c e N rather than on its sign, hence it remains 
almost unaffected by disorder. However, the oscillator 
wavefunction which is given by J2n e c n \N)ph> depends 
also on the sign of c e N and is affected by disorder. 

In Fig. 5 we plot | e (l|G)| 2 and | e (2|G)| 2 , which give 
the probabilities that the polaronic charge carrier in the 
ground state (|G)) lies at site 1 and site 2, respectively. 
For the ordered case these probabilities are same (=0.5) 
for any value of the coupling. With increasing disorder 
strength these probabilities differ from 0.5; the site with 
lower site-energy has higher probability of being occu- 
pied than the other. This difference in the occupancy in- 
creases rapidly during polaron crossover as a consequence 
of rapid suppression of the polaronic hopping (t e ) with in- 
creasing g + in the crossover region. For strong-coupling 
regime | e (l|G)| 2 becomes almost 1 while | e (2|G)| 2 ap- 
proaches to zero showing localization of the polaron at 
site 1. This localization as well as the correlation function 
X are shown in Fig. 5(a) and 5(b) for e<j = 0.5, 0.1 and 
0.001 to show that the polaronic crossover precedes the 
localization for t/w = 0.5. In Fig. 5(c) similar plots are 
presented for different values of t/w < 1.5 and td = 0.5. 
With increasing t the polaron crossover and the localiza- 
tion occur at higher values of e-ph coupling as expected. 
However, the localization is abrupt but continuous for 
larger values of t while the polaron crossover is still very 
smooth for t/w <1.5. 

It may be mentioned that when t/w >2 and td >0.5 we 
find that the polaron crossover and the localization occur 
almost simultaneously. However, for such large values oft 
our perturbation method does not converge properly in a 
small range around intermediate coupling where polaron 
crossover takes place. 

In Fig. 6 we plot x an d | e ( 1 1 G) | 2 as a function of 
disorder strength for different values of e-ph coupling 
with t/w=0.5. For very weak coupling (g + = 0.1) the 
polaron crossover is induced by disorder potential and 
is very smooth. Here the localization and the polaron 
crossover occur almost simultaneously but none of them 
is complete even at = 2. For g + = 0.4, which is also 
in the weak coupling range, the localization follows the 
polaron crossover and a large value of disorder strength 
is required for localization. For intermediate coupling 
(<7+ = 1) the localization takes place in the small polaron 
region and 95% of the localization is achieved within 
td/t ~ 1. For higher values of coupling the polaron would 
be a small polaron even in absence of disorder and local- 
ization would occur for a small value of disorder strength. 
Our study points out that when the e-ph coupling is in 
the intermediate range, localization may be achieved with 
a disorder strength of the order of half (electronic) band- 
width of the system; this has also been pointed out in 
Ref.[2j| in the context of manganites. For small e-ph 
coupling, disorder strength larger than the bandwidth is 
required for localization. 

IV. Two-site Holstein model with double 
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FIG. 6: Plots of x an d U(1|G)| 2 versus td for different values 
of g+ with t/u) = 0.5. \ is denoted by different symbols and 
| e (l|G}| 2 by broken line. 

exchange: effect of disorder 

The relevant Hamiltonian for studying a two-site dou- 
ble exchange Holstein model in presence of antiferro mag - 
netic interaction between core spins is given by 0, l24j 

H = e^ii + e 2 n 2 cos (9/2) (c f 1(J c 2a + cl a c la ) 

a 

+ gLuJ2 n *4h + bl) +u;J2 b l b * + ( 22 ) 

where Si, S 2 represent the local core spins (for man- 
ganites it is the spin of t 2g electrons) at sites 1 and 2, 
respectively and 9 is the angle between them. J is the 
superexchange antiferromagnetic interaction between the 
neighbouring core-spins S. The transfer hopping inte- 
gral (t) of the itinerant electron is modified to icos(|) 
because of the double exchange process which originates 
from strong Hund's coupling between the spins of the 
core electrons and itinerant electron |25j | . Here we would 
treat the core spins classically. For manganites the core 
spins have S = 3/2. However, for small values of J/t, 
the qualitative behavior of the phase diagram of the two- 
site Holstein-double exchange model does not depend on 
the value of the spin or hopping as observed in Ref . • 
Furthermore, D. M. Edwards and his group [23J pointed 
out that for such models the resistivity and transition 
temperature T c do not vary much with S, so that clas- 
sical spin is a convenient approximation to S=3/2 spins. 
Considering the out-of phase phonon mode which only 
couples with the electronic degrees of freedom and treat- 
ing the spin classically, we obtain the MLF transformed 
Hamiltonian [24| as 

n 

H c i = Lud^d + (ei — e p )ni + (e 2 — £ P )n 2 — tcos(-) 

c\c 2 exp(2A(d t - d)) + 4 Ci exp(— 2A(d^ — d)) 



+ uj(g+- \)(n 1 -n 2 )(d + d^) + JS 2 cos9 (23) 

In our previous work |24| we studied the above Hamilto- 
nian for a single polaron as a function of e-ph coupling for 
the ordered case (£1=62) using perturbation theory with 
the variational MLF basis. We found that the nature 
of the ferromagnetic (FM) to antiferromagnetic (AFM) 
transition as well as that of the polaronic state depends 
on the relative values of J and t. For small values of 
JS 2 /t the magnetic transition does not coincide with the 
polaronic crossover and a FM small polaronic state exists 
between a large polaronic FM state and extremely small 
polaronic AFM state. Similar phase diagrams have also 
been observed for both adiabatic and antiadiabatic limits 
in Ref. [lj| for small values of JS 2 /t. Here we will ex- 
amine the effect of site-energy disorder on such polaronic 
state and crossover. We have followed the procedure of 
our previous work p3 | to find out the ground state prop- 
erties of the double exchange Holstein model with site 
disorder. 

In manganites the ratio of the site energy disorder po- 
tential to the half band-width, when calcium is doped 
in LaMnC"3, is about 0.4 [23J. The value of the disorder 
strength used in the present work is of that order. It 
may be mentioned that previous studies |2(j| with a dou- 
ble exchange model have shown that the metal-insulator 
(MI) transition due to off-diagonal disorder (associated 
with random spins in the paramagnetic phase) requires 
also a large diagonal disorder strength. This led to the 
conclusion that disorder cannot alone account for the MI 
transition in manganites and an e-ph coupling of inter- 
mediate range is needed for this purpose [23| . 

In Fig. 7(a) and 7(b) we plot the angle (9) between 
the core spins, the correlation function x = { n i{ u i ~ 
u 2 )) /2g + (ni) and the kinetic energy {Ek.e) as a function 
of g + to show the FM-AFM transition, polaron crossover 
and the polaron-delocalization energy for both the or- 
dered and disordered cases. In presence of disorder, the 
FM-AFM magnetic transition and polaron crossover oc- 
cur at lower values of g+ . The reason behind this behav- 
ior is that disorder effectively reduces the hopping, hence 
favors formation of small polaron and the AFM phase. 
In the AFM phase, the polarons are small polarons with 
almost vanishing kinetic energy. However, the FM phase 
may have a large polaron character or a small polaron 
character depending on values of g+ . For low values of t 
(=0.7), it is difficult to distinguish between the regions of 
FM large polaron and FM small polaron since the change 
in x is small and the curves for both x and Ek.e are very 
smooth for disordered case. In Fig. 7(b) we consider a 
higher value of t (=1.5) with JS 2 /t = 0.05 as relevant 
for manganite systems |27j. Here a crossover from the 
FM large polaronic state to a FM small polaronic state 
with reduced kinetic energy is seen for the disordered 
case, before the transition to the AFM state. So the 
site disorder does not make any drastic change to the 
qualitative features of the ground-state properties of the 
double-exchange Holstein model but smoothens the po- 
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FIG. 7: Variation of 0, \ = ~ < m(«l — 1*2) > /2g+(n\) and 
Ba's with g+ for = and 0.5 in absence of magnetic field 
for (a) t = 0.7 and JS 2 = 0.05, and for (b) t = 1.5 and JS 2 = 
0.075. (c) The magnetoresistance MR=(E K E(h)-EKE(0))/h 
as a function of g+ for t = 1.5, JS* 2 = 0.075 and ft=0.04, ft 
represents the applied magnetic field. 



laron crossover. In Fig. 7(c) we plot the change in Eke 
due to the magnetic field (h) as a function of g + . This 
quantity may be related to the magnetoresistance for a 
system in the thermodynamic limit as pointed out in Ref. 
|24j . In general, the Eke is a measure of delocalization 
of the polarons and its change with the magnetic field 



gives field-induced delocalization of the polaronic charge 
carriers. In Ref. [24[ we reported for the ordered case 
that for JS 2 /t =0.05 the change in Eke due to the field 
has a broad peak around the FM-AFM transition. Wc 
find that disorder makes the peak more broader but the 
value of the change in Eke remains almost the same. 
We believe that the magneto-resistance of similar model 
system in the thermodynamic limit would show similar 
qualitative features. We do not present here the results 
for larger values of JS 2 /t where the magnetic transition 
coincides with the polaron crossover and the change in 
the Eke shows a sharp peak at the FM-AFM transition 
[24| as we find that the site disorder does not change that 
behavior. 



V. Conclusions 



To summarize, we have presented the results on the 
two-site single polaron Holstein model in presence of a 
site energy disorder e<2 which appears as a difference in 
site energies. With increasing the retardation between 
the electron and associated deformation becomes weaker 
and the polaron crossover occurs at lower values of e- 
ph coupling. The polaronic kinetic energy is suppressed 
appreciably with disorder in the range from weak to in- 
termediate couplings. But in the strong-coupling region 
where the hopping through multi-phonon process is dom- 
inant, kinetic energy is almost independent of td- For the 
oscillator wavefunction, a broad peak in strong-coupling 
region is observed for the ordered case (ed—0) in addition 
to the main peak. In presence of disorder (e^ =1.0), this 
feature disappears. We find that even a small disorder 
(ed=0.001) can localize the electron in presence of e-ph 
coupling. The polaron crossover precedes the localiza- 
tion in the non-adiabatic regime for low to intermediate 
values of disorder (ed < t) . 

For the double exchange Holstein model both the mag- 
netic transition and polaron crossover shift toward lower 
values of e-ph coupling with increasing disorder strength. 
The qualitative features of the ground-state phase dia- 
gram does not change much with disorder. The magneto- 
resistance for the model shows up in a broader region of 
parameter space in presence of disorder. 
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